The Ising M-p-spin mean-field model for the structural glass: 
continuous vs. discontinuous transition 



ii 



F. Caltagirone 1 , U. Ferrari 1 , L. Leuzzi 1,2 , G. Parisi 1,2 ' 3 and T. Rizzo 11 
Dip. Fisica, Universtta "Sapienza", Piazzale A. Moro 2, 1-00185, Rome, Italy 



2 IPCF-CNR, UOS Rome, Universtta "Sapienza", Piazzale A. Moro 2, 1-00185, Rome, Italy 
3 INFN, Piazzale A. Moro 2, 00185, Rome, Italy 
(Dated: October 26, 2010) 

The critical behavior ol a family of fully connected mean-field models with quenched disorder, the 
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is straightforwardly extendable to finite dimensions in any geometry. 
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I. INTRODUCTION 

Since the work of Kirkpatrick, Thirumalai and 
Wolynes 1-6 a certain set of mean-field spin-glass 
models have been shown to own the salient proper- 
ties of the behavior of structural glasses. In particu- 
lar, these models display dynamic equations that are 
equivalent to those predicted by the Mode Coupling 
Theory (MCT) 7 ~ 9 above the so-called mode coupling 
temperature T mc where ergodicity breaking occurs in 
that theory. Moreover, two kinds of transition are 
predicted: a dynamic one at = T mc and a ther- 
modynamic phase transition at a lower T, often re- 
ferred to as Kauzmann transition. Mean-field mod- 
els exhibiting structural glass features are charac- 
terized by multi-body microscopic interactions and 
their thermodynamics is self-consistently described 
by implementing a discontinuous Replica Symmetry 
Breaking (RSB) Ansatz (usually one step: 1RSB). 

The dynamic transition is due to the presence of 
a large number of metastable excited glassy states, 
represented as local minima in the free energy land- 
scape in the configurations space, growing exponen- 
tially with the size N of the system. In the mean- 
field approximation, barriers between minima grow 
with the size, so that, in the thermodynamic limit, 
the relaxing dynamics to equilibrium of the system 
at T < Td remains stuck forever inside the first 
"meta" -stable state where it ends up in. In real 
glassy systems, however, there is a slow dynamics 
occurring through activated processes and this dy- 
namic arrest is an artefact due to the mean-field 
approximation. In finite dimensions the glass tran- 
sition occurs because at some glass temperaure T g 
the time-scales of observation are shorter than the 
characteristic time-scales of the slowest structural 
a processes taking place in the glass-former sam- 



ple. Metastable states really have a finite time- 
life, even though (much) longer than the experi- 
mental time of observation. The effect of activated 
processes starting from spin-glass 1RSB mean-field 
models has been analyzed, e.g., by working at fi- 
nite N in the fully connected Random Orthogonal 
Model (ROM) 10-13 and finding a glass behavior, sim- 
ilar to the one observed in computer glasses, cf., e.g., 
Refs. 14 - 15 . 

As we mentioned, another property occurring in 
the glass-like mean-field models (see also Refs. 16-18 ) 
is a thermodynamic transition between the super- 
cooled liquid (below Td) and a thermodynamically 
stable glass. This occurs with a jump in the order 
parameter, but without discontinuity in the internal 
energy (no latent heat is exchanged). This mixture 
of first order and continuous phase transition in pres- 
ence of disorder has been termed random first order 
(RFOT). 6 

One of the most accredited theories, the Adam- 
Gibbs-Di Marzio entropic theory 19,20 predicts the 
existence of a thermodynamic transition to an ideal 
glass phase, the so-called Kauzmann transition. The 
Kauzmann temperature is generally associated with 
the asymptote of the Vogel-Fulcher law 21,22 of the 
relaxation time and, thus, related to the transition 
one might have in a infinitely slow cooling of a never- 
crystalizing glass-former. Because of the impossibil- 
ity of experimental measurements of glass relaxation 
in those conditions, the very existence of the Kauz- 
mann point and the nature of that transition is still 
a matter of debate. 

Attempts to follow the properties envisaged in 
mean-field models in realistic systems have faced 
the problem of finding a proper way of embedding 
the model microscopic features into a given finite 
dimensional geometry (e.g., 3D cubic lattice) with- 
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out altering the discontinuous nature of the transi- 
tion. So that one can actually falsify the hyphothesis 
of RFOT in finite dimensional systems. Indeed, in 
Ref. [23] a generalization of the p = 3-spin model 
with M — 2 Ising spins on each site was numerically 
studied on a D — 4 hypercubic lattice finding evi- 
dence for a continuous phase transition. The same 
continuous behavior was recently found, already in 
the mean- field regime, in the same p = 3, M = 2 
model in a D = 1 chain on a "Levy lattice". 24 

Starting from this observation, that the RFOT 
becomes continuous in finite D, the work of Moore, 
Drossel and Yeo 25 ~ 27 shows that this is equivalent 
to the critical behavior of the Edwards-Anderson 
model in a field, where the transition line is called de 
Almeida-Thouless (dAT) line. Applying droplet the- 
ory (that rules out the existence of a dAT line out- 
side the limit of validity of mean-field theory) it is, 
thus, inferred that no thermodynamic random first 
order transition can occur in real structural glasses. 
The issue of the existence of a dAT line in finite di- 
mensional spin-glasses will not be addressed here. 
For recent bibliography on that subject see Refs. 
[28-33] and references therein. 

In the present work we will focus on deriving 
a mean-field class of models, to whom Ising p- 
spins 16 ' 34 belong, whose critical behavior shifts from 
continuous to discontinuous in a controlled way. The 
aim is to clarify why the finite dimension extensions 
of mean-field glasses studied so far do not display 
RFOT and to devise mean-field models whose dis- 
continuous critical nature can be conserved also be- 
yond the limit of validity of the mean-field approx- 
imation. The model consists of N sites, each one 
containing M spins interacting with spins on other 
sites in p-uples. We will see how, changing p and 
the number M of spins living on a single site, it 
is possible to move from systems displaying a sec- 
ond order phase transition to systems displaying a 
random first order transition, that is, yielding both 
a dynamic and a Kauzmann-like phase transition. 
The finite dimension counterpart of the model under 
probe can be easily achieved since the p-spin inter- 
action is always exchanged between two sites, e.g., 
nearest neighbors on a d-dimcnsional (hyper)cubic 
lattice. 

We mention that moving from mean-field to fi- 
nite dimensions, also standard Ising p-spin and Potts 
models might conserve the random first order nature 
of the transition and keep reproducing basic features 
of structural glasses. Even though it is not straight- 
forward to conceive a short-range finite dimensional 
Ising p-spin, the Potts model can be easily defined 



on a hyper-cubic lattice. Nevertheless, no numerical 
evidence has been collected so far for a discontinu- 
ous RFOT in disordered Potts models with number 
of states ppotts = 5,6, 10 35 ~ 37 and, actually, we found 
no argument to infer that in the finite dimensional 
lattice case the pp tts — > oo limit can be kept under 
control. 

On the other hand, the model considered in the 
present work has the advantage to reduce to an ex- 
act mean-field model for the RFOT as M — > oo even 
in finite dimension (and finite size), for any values of 
p. Moreover, we can work out a sufficient criterion 
to determine the smallest value of M above which 
continuous transitions cannot occur. 

The manuscript is organized as follows: in Sec. II 
we will study the statistical mechanics of the model; 
in Sec. Ill we show that the large M limit cor- 
responds to standard p-spin and in Sec. IV, ex- 
panding near criticality, we build the correspond- 
ing field-theory, compute the coupling constants and 
study the relevance of terms competing for continu- 
ous/discontinuous transition. In Sec. V we present 
our conclusions. 



II. THE MODEL 

The model consists on N sites, each one hosting a 
set of M spins. Two sites interact through a p-body 
interaction involving spins belonging to the two sets 
of M spins. The Hamiltonian reads 



«=-e e 



(i) 



M69 



where < x, y > indicates the sum over all couples of 
sites and g(x, y) are all the possible p-uplets among 
the 2M spins, with an exception if p < M: those 
p-uplets completely pertaining to a single site are 
excluded. This choice actually defines our model 
when p < M, as we will discuss in the following. 

The disordered interactions are Gaussian i.i.d. 
variables, with distribution: 



P(J a ) 



(2) 



where, to provide the right thermodynamic conver- 
gence of the free energy, the variance scales like 
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A. Free energy and order parameters 

Replicating n times the system we compute the 
average over quenched disorder of the replicated par- 
tition function: 

- l,N 
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Explicitly separating those spins belonging to site 
x from those on site y one can obtain a general 
expression for the partition function valid both for 
p > M and p < M: 



Z n = Tr 



{s(x),s(v)} eX P 
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For p > M the sum over fc runs from p - M to M; 
in the case p < M the sum over k runs from 1 to 
p— 1. In principle, it might be possible to include an 
extra term due to self-interaction: p out of M spins 
interact on a single site ("a single site standard p- 
spin"). As already mentioned, in the present work 
we will consider a model without site self-interaction. 
We now introduce a set of multi-overlaps between k 
spins on the same site x in two replicas: 
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with 



fc = p - M, . . . , M if p > M 
k = 1, ... ,p- 1 if p < M 



By means of multi-overlaps we can write the repli- 
cated partition function Eq. (7) as 
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where the parameter C, proportional to minus the 
paramagnetic free energy, reads 
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Introducing the integral representation for the 
delta functions in Eq. (8) one obtains: 
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The stationarity equations in A and Q are 
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Substituting the saddle point value for A in the ef- 
fective action we obtain 
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The physical meaning of the overlap matrix at 
saddle point value is the usual one and, more pre- 
cisely 
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III. LARGE M LIMIT: STANDARD p-SPIN 



we obtain, for the free energy Eq. (16) 
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The saddle point self-consistency equation w.r.t. 

Q(p-k) yields 
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Substituting Eq. (21) in Eq. (20), we obtain the 
expression 

G(q,X)/M = --J2J2 kUp-kY q " b + Xabqab 
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that is, the standard formal free energy of the fully 
connected Ising p-spin model: 
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with 



For large M, neglecting diagonal terms in the sum 
over . . ,ik, in Eq. (16), the logTr term can be 
rewritten as 
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Performing the saddle point for large M, rather 
than N, and introducing the auxiliary parameter 
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IV. ANALYSIS OF THE CRITICAL POINT 

Our aim is to find the transition point and to 
study its thermodynamic nature as M and p are 
changed. In particular, we will verify that, at given 
p (vice-versa M) there are threshold values of M 
(resp. p) beyond which the transition switches from 
continuous to discontinuous. 
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First, to identify the critical point we expand the 
stationarity equation (14) to first order in Q^ b \ ob- 
taining: 
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There are "multi"- critical temperatures for the 
"multi" - overlaps, whose expressions read 
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The largest critical temperature is obtained for k = 
p/2 if p is even, and for k = (p + l)/2 , (p — l)/2 if 
p is odd. The overlap corresponding to the smallest 
(i c (slightly above /3 C ) is non-zero and of order t oc 
(T c — T)/T c , while the others are at least of order 

T 2 . 

Proceeding to the second order expansion of Eq. 
(14) we have 
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we will focus only on the equations for the overlaps 
corresponding to the largest critical temperature, cf. 
Eq. (27), that is, on the terms of the type 



Q (3 fW 2) 

^ab "» ab 



for even p 



In Eq. (28) each spin in each replica has to be 
matched with another one in another replica in or- 
der to get a non-zero result from the trace. At sec- 
ond order we are, thus, left with only two kinds of 
possible matching, yielding terms: 



We will see how, depending on the parity of p even 
the multiplicity of such terms will change, leading to 
different expressions of their coefficients as functions 
of p and M. 

Using the above results, Eq. (16), approximated 
to the second order in Q, can then be written as 



a, b 
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where Q ab stays for Q y ab . 

As already noticed by Gross, Kanter and 
Sompolinsky 17 in the Potts model (threshold was 
PPotts = 4 colors) and in Ref. [38], it can be shown 
(see App. A) that if the ratio between coupling con- 
stants on the nonlinear terms is larger than one the 
phase transition cannot be continuous. We will now 
proceed to the computation of the coupling con- 
stants for the M-p Ising spin model. Since, as al- 
ready mentioned, the computation of the third order 
coefficients will yield different functional expressions 
depending on the values of p, we have to distinguish 
between four cases: 
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and we will analyze them separately. 
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More specifically, we are interested in the terms of 
the series at the r.h.s of Eq. (28) with k = I = m = 
p/2, Hp is even, or with k,l,m = if p is odd. 

It is interesting to notice that we would have 
the same physics considering a model in which p/2- 
uples on each site interact with p/2 on another site 
(p even) or (p + l)/2-uples on a site interact with 
(p— l)/2 on another site (p odd). 



A. Even p and p/2, p — 4a 

The only surviving term in the sum over I and m 
in the r.h.s. of Eq. (28) is for I = m = k = p/2. The 
trace term turns out to be 
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The squared term is 
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B. Even p, odd p/2, p = 4a - 2 



odd. For p = 4a — 1 we obtain 
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In the r.h.s. of Eq. (29) only the coefficient in 
front of the Tr Q^ b term survives, whereas w 2 = 
always. The ratio is 



w 2 



= 



(34) 



According to the small Q expansion, Eqs. (28), (29), 
when p is even and p/2 is odd the transition at the 
largest critical temperature l//3 c (p/2), cf. Eq. (27), 
turns out can to be continuous, no matter how many 
spins M stay on each site. This might appear at con- 
trast with the large M limit leading to Eq. (23) that 
is equivalent to an Ising p-spin mean-field model for 
any p > 2. There, however, no pcrturbative expan- 
sion was carried out, while Eq. (34) is the outcome 
of an expansion for small overlap values that can- 
not help in identifying discontinuous transitions with 
large jumps in Q. Indeed, the condition expressed by 
Eq. (A7) is sufficient but not necessary to rule out 
a continuous transition. 



C. Odd p, even (p+l)/2, p = 4a - 1 

When p is odd we have to deal with two relevant 
overlaps Q^-^ - ^ and Q^ 1 ^^ and one critical temper- 
ature. In order to determine the coupling constants 
of the cubic terms one, thus, has to diagonalize a 
2x2 matrix. In App. B we report the details of the 
computation leading to: 
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4Mp- 3 / 2 vpTT V 



for the coefficient of the cubic trace term in the ac- 
tion, Eq. (29). The expression for the coupling con- 
stant depends, further, on (p + l)/2 being even or 



D. Odd p, even (p - l)/2, p = 4a - 3 

In this last case the coupling of the trace cubic 
term is still given by Eq. (37) and the second non- 
linear coupling constant is expressed as: 



w 2 (p,M) = 



yielding the ratio: 
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We now have a complete description of the critical 
behavior of the M-p system. Already at the mean- 
field level, in order to have a discontinuous transition 
a p > 2 interaction between spins in not enough. 
We find that for each given p one needs a minimal 
number of spins Md; sc on each site in order to have 
a random first order phase transition, corresponding 
to the lowest integer M for which w 2 /w\ > 1; cf. 
Eq. (33), (37), or (39) depending on the parity of p 
and (p+ l)/2. 

In table Tab. I we report some values of the ratios 
for systems with small p and M. In Fig. 1 we plot 
the Mdi S c(p) behavior. 



V. CONCLUSIONS 

In the present work we have performed an analytic 
computation of the critical behavior of a mean-field 
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TABLE I: Ratio values for small p and M around the 
threshold 1 
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FIG. 1: Lowest integer values of M at given p, for which a 
discontinuous transition is certainly expected (sufficient 
condition to have RFOT is M > M d i sc ). 



p-spin model that can display both a random first 
order and a continuous phase transition tuning the 
number of spins present on each site. The effective 
action is represented by Eq. (29) where the first two 
relevant terms are third order in the overlap param- 
eter Q. The ratio of their coupling constants has 
been computed, thus providing the threshold val- 
ues of Mdi S c (p) separating models with continuous 
from models with discontinuous critical behavior. If 
M > Mdisc this is sufficient to guarantee the discon- 
tinuity of the transition, though it is not necessary. 
A RFOT can also occur at lower values of M that 
cannot be identified with our probe based on the 
perturbative expansion for small Q, cf. Eq. (29). 

The particular case studied in Ref. [23], (M = 2, 
p = 3), yielded numerical evidence for a continuous 
phase transition in dimension four. This is consis- 
tent with the value of the w^/wx — 3[1 — 1/V2] = 
0.87868 as computed in the mean-field theory, cf. 
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Tab. I and Ref. [39]. The same applies to the 
model recently studied in Ref. [24], a one dimen- 
sional (M,p) — (2, 3) model on a Levy lattice 40 . 

The continuous-discontinuous cross-over is like the 
one found in Potts 17 varying the number of colors, 
and in the spherical p-spin varying an external mag- 
netic field 18 . The advantage of the present model is 
that it can be easily represented in finite dimensions 
on lattices of given geometry, e.g., on a cubic lattice 
with short-range interactions, and that it always dis- 
plays a RFOT in the M — > oo limit for all lattices, 
both finite dimensional and mean-field like. 
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Appendix A: Threshold value for W2/1C1 

Starting from the self-consistency equation for 
small Q's (implying a continuous transition in Q), 
cf. Eq. (29) where also the quartic term is consid- 
ered, 

TQab + w 1 {Q 2 ) ab + w 2 Q 2 ab + yQl b = (Al) 
we have, in the RSB Ansatz, 



rq(x) — W\ 2q{x) I q(s)ds 
'0 



(A2) 



+ / [q(x) - q(s)} 2 ds + w 2 q(x) 2 + yq{xf = 



(A3) 



Deriving once Eq. (A2) w.r.t. x one has 



q'(x) \ T ~ 2wi 



q(s) ds + xq{x) 



+2w 2 q(x) +3yq{x) 2 *j = 

If q'{x) ^ 0, deriving a second time w.r.t. x, one 
finds 

q' (x) [-W!X + w 2 + 3yq{x)} =0 (A4) 

If y > 0, then the overlap function around criticality 
can be written as: 



x < x\ 

q{x) — ^ j^(wix — u> 2 ) X\ < x < X2 (A5) 
q(l) x > x 2 
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where, for continuity, x\ = Wijw\ and 
t 2 



9(1) 
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p = ui!(l — x x ) (A6) 



As a consequence, in order to have a continuous 
transition it must be 



W2 
Wi 



< 1. 
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The argument for the threshold value of x still 
works also if y < in Eq. (Al). In that case, rather 
than a continuous function, we simply have a 1RSB 
step function for q(x) = 9(x — xi)q, with 
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Appendix B: Coupling constants with odd p 

When p is odd we have to deal with two relevant 
overlaps and one critical temperature. The second 
order equation, Eq. (28), has the structure: 



AQ ab = F({Q}) 



(Bl) 



where Q a6 = {Q^ b M \--- ,Q^}. Diagonalizing 
A -)• V A = V~^AV one obtains: 

V^AVV^Qab = V^FdQat,}) (B2) 

Introducing new variables & a b, linear combinations 
of Qab, the above expression can be rewritten as 



V A @ ab = V-'FdV&ab}) 



(B3) 



Rearranging the entries in a proper way, A can be 
written as a block matrix of 2 x 2 elements per block, 
and each block can be diagonalized separately, with 
eigenvalues 



A (k±) = 1 ± 02 ^f( k )f(p-k) (B4) 
and eigenvectors: 



,k± 



T; 1 



Wf(k) J 



(B5) 



For each block of the matrix, labeled by k, the eigen- 
vector matrix and its inverse, thus, are 



V = 



■P = 



2y/f(p-k) 2y/f(p-k) 
1 1 

2v7(fe) 2v7W 



V /(p-fc) -v/JW 



with 



and 



/(*0 



1 



M 2fe_1 V fe 



(B6) 



(B7) 



(B8) 



e (fe+) - V/(p-*)Q (fc) - v7WQ (p_fc) (B9) 

© (fe_) = y/f(P~k)Q {k) + \ffik)Q [p ~ k) (BIO) 
In the present case, since p is odd, the only overlaps 
we need to consider are Q( £ i i ) and QW. Their 
self-consistency equation can be written in the form: 



A 



Vab . 



with 



.4 



-/3 2 /(^) 



^z_i(Q) 

F £±i (Q) 



-/3 2 /(^) 
1 



(Bll) 



(B12) 



The functions i r ( p ±i)/2 are two polinomials of degree 
two in all the Q^''s. However, as mentioned above, 
in order to study the nature of the critical behav- 
ior (continuous or discontinuous) we only need the 
terms relevant at the highest critical temperature, 
cf. Eq. (27), and we, thus, set to zero all the over- 
lap matrices except for Q^^a - ) and Q*- 2 ^ - -*. 

Depending on the parity of {p + l)/2 the relevant 
terms contributing to the nonlinear couplings in the 
action Eq. (29) differ. We now consider the two 
cases separately. 



a. (p + l)/2 even 

If p = 4a— 1 with a e N the functions on the r.h.s. 
of Eq. (Bll) read: 



cb 



8* /M\ A ( E±i) 

\ z / C—l 

(B13) 



Applying the diagonalization transformation de- 
scribed above, cf. Eqs. (B1)-(B3), one finds 



^ 2 ) (C^) (B19) 

(B14) 



(B20) 

where the eigenvalues, cf. Eq. (B4), are 



b. (p — l)/2 even 



A (±) = 1 ± /?V/((p-1)/2)/((p+1)/2) 



Since the F's depend on the Q's, we have to apply 

the inverse transformation to get equations in terms 
If, otherwise, p = 4a + 1 with a e N, one obtains of the Q , s _ The eigenvalue A (+) is always pos i t i vej 

so that 8^ + ) plays the same role of the "non criti- 
p p _ 1 — ft ( M \ q^T^q' 2 ^ ca l" overlaps and can be put to zero. The inverse 

Miv~i \ ^tt- I ^-^ "" transformation, thus, reduces to 



m\(^\(m-^\ ( ^ ^ 9 (-) 

P-1 ) \ p-1 ) \ p-1 j y^ab ) Wab 



4 / \ 4 



+ 2M fp-|(E^i)(4i)( P+3 2 )( Q ^ 2) ) 



^ ( M\(*±\(M-*i\ r „r-^)\ 2 n( - } ( B21 ) 



(B15) V^ 2 ^) 
so that Eq. (B20) decouples in 



'cb 



, g fM N /^N /M - ^ W>Q^> +W2(P ' M) (9 - ))2 ' (B22) 

Af § \2±1 y ^£^1 y ^ 2^1 J Wab ^ ab The constants w\ and w 2 depend on p and M. The 



(B16) 



Computation of the coupling constants (B23) 

The formula for w 2 changes depending on the parity 
In order to decouple Eqs. (B11)-(B12) we specify of (P + !)/ 2 - For cvcn (? + 
the two new variables, Eqs. (B9)-(B10), for k = < 

{p - 1)/2: MP,M) = ^-^^( £ _ i ) (B24) 



(B17) 



(id-*? 



p+i i i p+i 

4 / \ 4 

E±l\ /M-2±l 

2+1 J 1 p+1 

4 / \ 4 



e«;> - J, (£±i eiP + . / ^ elf 2M - P+1 ^ /m - «±i 



+ 



(B18) 1+P V^/V 2? 
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If (p + l)/2 is odd it reads: 
1 



w 2 (p, M) 



M 

MI'-* \/ KH 1 



(B25) 



2 

P-l 
4 



4 



+ - 



1+P 



2M-p + 1 V 2 ^/ V 2j r- 



+2 



£±1\ /M- 2±1 



2 

P-l 
4 



P-l 
4 
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